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O ■ Abstract 

Q> ' A model of quark confinement based on a singular solution of classical YM 

r~| ' equation is proposed. Within the framework of this model we have calculated 

Qh! hadron masses that correspond to ground state configurations of quarks. Our 

Q^' results are in agreement with the experiment data with accuracy 3-7 percents 

^ ■ for all hadronic masses except those of light pseudoscalar mesons. 
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1 Introduction. 

Bag model is one of the first attempts to understand physics of hadrons in terms 
of quarks. In its simplest form it was formulated in the end of 60 's and became 
well-known after the works of MIT group in the middle of 70's. So far many mod- 
ifications of bag model were proposed 0. It appeared that bag model gave rather 
satisfactory description of masses, magnetic moments and many other properties of 
hadrons . However, in spite of these successes, bag model never was a very popular 
one. The origin is clear. In contrast to Standard model pretending for description 
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of hadron physics by proceeding from some fundamental principles, bag model is 
based on a pure phenomenological assumption that quarks are confined inside a cer- 
tain sphere. Up to now it is absolutely unclear how to derive such an assumption 
from QCD and just so nowadays bag model is usually considered as a rather crude 
phenomenological model, and nothing more. 

In the present paper we will try to show that , nevertheless, there exists a way to 
derive a kind of bag model from QCD. 

Our basic assumption is that quark in zero approximation moves in a certain 
effective YM potential, that is a solution to classical YM equations with singularity 
on the sphere. Such solutions were discovered in 70's in papers 0. More recently 
singular solution of YM equations was found as an analogue to Schwartzchild solution 
in general relativity . In the latter paper a remarkable analogy between YM theory 
and general relativity established in the works was used. (See for further 
references). Later analogous solutions for Yang-Mills-Higgs equations were found 0]. 
Solutions with singularity on the torus and cylinder were investigated in ||1CI|| . 

Our model obviously can be derived from QCD by quantization in neighborhood 
of such a singular classical solution as zero approximation. Further corrections, in 
principal, can be also obtained in a systematic way. But in present paper we restrict 
ourselves to investigation of zero approximation. Namely, we will investigate the 
motion of quark in such field and evaluate mass spectrum of ground hadron states. We 
will show that our model gives quite satisfactory description of all ground hadron state 
except those of light pseudoscolar mesons. The discrepancy between our results and 
experiment data for masses of light pseudoscolar mesons is not a surprise. Indeed, it is 
known that many features of light mesons physics are related to spontaneous breaking 
of chiral symmetry. But the latter phenomenon is not taken into consideration within 
framework of our bag model at all. So one cannot hope to describe light meson mass 
spectrum within framework of the present simple version of our model. This problem 
needs further investigations. 

The paper is organized as follows. In section 2 we describe spherically symmetric 
solutions of YM equations with singularity on the sphere. In section 3 we investigate 
motion of Dirac particle in such field. In section 4 we calculate the mass spectrum of 
ground hadron states by quantization in the neighborhood of this classical solution. 
In the last section we discuss results obtained. 



2 Classical solutions of YM equations with singu- 
larity on the sphere. 

Let us consider four dimensional SU(2) YM theory. Substituting the well-known 
Wu-Yang ansatz [|11[] 

A^, = 0, A1 = ea^~{l - H) (1) 
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in YM equations, one get 

r^H" = H{H^ - 1) (2) 

It can be proved that there are only two types of solutions of equation (2) that 
are singular on some sphere and regular at r = (see Fig.l). 
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Figure 1: The singular solutions of Wu-Yang equation (2), i? = 1. 
The asymptotics of the first one are: 

R-O; Hi{r)r^l + Ci{^f,r-.0 (3) 



R-0; H2{r) ^ -1 + C2{^)\ r ^ (4) 

where i? > is an arbitrary constant, Ci ~ 3.038 and C2 — 9.448. 

It appears that satisfactory agreement between experimental data and results 
obtained from our model can be achieved only for solution with asymptotics (4). 



HAr) 

^ ^ 1 -r/i? 

and of the second one are: 
1 — r/R 
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So in what follows we will suppose that functions H{r) is defined by eqs.(2),(4). 

Starting from these singular solutions of SU{2) YM equations, one can construct two 
different sets of singular solutions of SU (3) YM equations that correspond to two 
non-equivalent embeddings of the algebra SU{2) in the algebra SU{Z). Namely, up 
to unitary equivalence, the first set of such solutions can be defined as 

A = \a^\^ + \a^\^ + \a^\^ (5) 

Zi 

and the second one as 

A = -A^\^ + -A^X'' + -A'-X' (6) 
2 2 2 ^ ' 

where A", a = 1, 2 • • • 8, are Gell-Mann matrices and A", a = 1, 2, 3 are defined by 
formula (1). 



3 Quark in Yang-Mills field with singularity on the 
sphere. 

Let us consider solutions of Dirac equation 

Q 

(i7°— + i7^Vj)* = (7) 

for the quark in potential (5) or (6). In the case (5) quark cannot be confined inside 
the sphere r = R because the third component of quark field satisfies free Dirac 
equation. So in what follows we restrict ourselves to consideration of the case (6). 

By virtue of spherical symmetry of the potential (6) Dirac Hamiltonian commutes 
with operators of angular momentum 

where (/") = (A^, A^, A'') are color isospin operators and are orbital angular mo- 
mentum ones. 

The solutions with definite total angular momentum and definite energy E can 
be represented as 




E+m 



(8) 
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where J^jm 

spherical spinors 



P,J+l/2 _ f~ V i( J+t/ ^-^+1/2 M-l/2\ f.J-1/2 _ / \/^Ty^-^-l/2M-l/2^ 

^ V 2(J+1) ^J+l/2M+l/2 / \^2(JTT)^"'-V2M+1/2/ 

and 

Vj = r- Xi— ni = — 

oxi or r 

One notes, that J = 1/2, 3/2, ■ ■ -. Substituting (8) in (7), one gets: 

Bt' = Y^{-(J - V2)(J + 3/2)^5+ + -Dt - (J - l/2)icr 
J + 1/2^ r r r 

-S{Cl - (J - l/2)fir)) (all J) 



Bi' = J^iiJ - 1/2)(J + 3/2)l5r - ^Di - (J + 3/2)ic+ 

-S{C^ + {J-l/2)Bt)) (J>l/2) 

1^1+' = -(J + l/2)-Dt + (J + 3/2)-B+ - (J - l/2)-Cr 

' ryi /yi 

+SD{ (all J) 

Dr' = (-f + 1/2) -/^r + iJ + 3/2) -c+ - (J - 1/2) -5r 

ry ry ly 

-SDt (all J) 

Cr' = 7:^((^ - 1/2)(J + 3/2)V - -Dt - (J + 3/2)^5+ 
J + 1/2^ r r r 



+^(fir + (-f + 3/2)Ci+)j (J>l/2) 

C+' = 7^(-(J - 1/2)(J + 3/2)ic+ + -Dr - (J - l/2)-5r 

+^(5+-(J-l/2)Cr)) (all J) 



(9) 



where £ = V E'^ — 'm'^, and exactly the same system of equations for functions i?^, 
C2, Df. So initial Dirac equations are separated into two identical sets of equations. 
This leads to degeneration of all energy levels. 

This is not surprises. Indeed, Dirac equations in external chromomagnetic field 
with property {A{—x) = —A{x)) are invariant under parity transformations . It is 
well-known ||12[ that this symmetry implies = 1 sypersymmetry of Dirac equations 
that, in turn, implies the doubling of energy levels. Thus all energy levels appear to 
be degenerate in parity. As a consequence, all hadrons must be degenerate in parity. 
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This is a puzzle of our model. In what follows, we postulate that quarks may be in 

states with only one parity. 

States with definite parity correspond to solutions of Dirac eq.(7) for which either 
Bf = = Df = or B2 = C2 = Df = 0- definiteness, we choose the second 
case. 

li J — 1/2 then only four functions -B/'(r), -Di"(r), C^(r) and -Df (r) survive in 
eqs.(9). Surprisingly, but four first order equations for the functions B^'^r), Df(r), 
Ci{r) and D^^r) reduce to one second order equation. Namely, if 

x+ = Dt + V2B+, x_^D+ - V2B+, 
y+ = + V2Ct, y- = - x/2C+ 
then simple analysis of eqs.(9) shows that either 

y_ = 

- ^ - ^)y^ = o ^^^^ 
x-^\{y+'-{i + ^)y^) 

or 

y+ = ^ 

y- +{£\-'-^ + ^)y- = (11) 

= liy-' - ii - ^)y-) 

Eqs.(lO) as well as eqs.(ll) has only one square locally integrable solution at 
r — R. It automatically satisfies boundary condition (12) 

y±{R) = (12) 

This means that quark cannot penetrate through the sphere r — R. Indeed, 
equations describing quark penetration through the potential barrier at r = i? must 
have at least two linear independent locally square integrable at r = i? solutions that 
correspond to two possible directions of quark motion (into the sphere and out of 
the sphere). Lowest eigenvalues £ = Sn/R of eqs.(10),(ll) with function H defined by 
eqs.(2),(4) are given in table 1. 

Energy levels can be can be express via numbers £n as 




En{R)^\^,+m^ (13) 



Let us consider the case J > 3/2. In this case the equations (9) have four linearly 
independent solutions that are locally square integrable at r = i?. Two of them vanish 
at r = i? whereas two other have asymptotics 
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7.800 


1.997 


10.920 


7.288 


14.074 


10.885 



Table 1: Eigenvalues S = Sn/R for eq.(lO) and (11). Data in the left and right 
columns correspond to eq.(lO) and (11) respectively. 
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(14) 



+ 0(1 - r/R) 



(15) 



/ BT \ 
BT 
Dt 
DI 

cr 
V ct J 

where Ci, C2 are constants. 

This means that quarks with J > 3/2 are not confined inside the sphere r — R. 
In particular, normal component of the carrent ji, — 'i/j^u'^ on the surface r — R is 
not vanished: 



\r=R 



C1C2 



In such situation the only possibility to confine quark with J > 3/2 sphere r 
is to impose ("by hands", as in usual bag models) the boundary condition 



\r=R 



= 



(16) 
= R 

(17) 



Obviously eq.(17) imphes C1C2 = and so either Ci = or C2 = 0. In other 
words, we must forbid either asymptotics (14) or asymptotics (15). We choose the 
condition Ci = 0. 

Results of numerical evolution of lovest eigenvalues S of the system (9) with 
boundary condition Ci = are presented in table 2. 
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J = 3/2 


J = 5/2 


4.165 


5.151 


6.977 


8.518 


7.949 


8.863 



Table 2: Eigenvalues S = Sn/R for eq.(9), . Data in the left and right columns 
correspond to J = 3/2 and J = 5/2 respectively. 

4 Hadron mass spectrum. 

Let us consider "partition function" 

Z = tr e 

where H is QCD Hamiltonian. "Partition function" can be represented as 

Z = J DAD^D'^e'-^o '^^^'^'^''^^^^''^^^^'^f'''-"'^^''^'^^^ (18) 

We assume that the main contribution in functional integral (9) is given by tra- 
jectories close to classical solution A^i of YM equations defined by eqs.(2),(4),and (6). 
Applying stationary phase method, one gets in zero approximation: 

Z= I e-^--(«)^n^e^[7°(n°^ + nV(^0-m,)]^P^C7 (19) 



zero 
modes 



where APBC means " anti-pereodic boundary conditions" and Eym{R) is classical 
energy of the field A^. The determinant in (19) can be easy evaluated in terms of 
positive eigenvalues Es{R,mg) (defined by eq.(13)) of Hamiltonian 

Z7°7V(A,0 - 7°m, 

that have been evaluated in section 3 (see, for instance, |T^). Substituting the result 
in (19), one obtains: 

Z = J e"*^^^*'*^'^^^^'!^''''"^"^^'"'''^^^ (20) 

1</Cs<ns zero 
modes 

where Hs = 4(2J + 1) is doubled multiplicity of the eigenvalue Eg (Remind, that 
multipl6icity of the eigenvalue Eg with given J is 2(2J + 1). Additional factor 2 
correspond to contribution of anti-quarks.). 
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By virtue of scale invariance of YM equations the measure of integration in (20) 
comprises integration with respect to R. So, applying again stationary phase method, 
one gets 



l<ks<ns 

where Rq is defined from an equation 



-^[EvMiR) + hE,{R, m,)] 



So hadronic masses are defined by the formule 



(23) 



R=Ro 



M'^'"' = Eym{Ro) + E E{Ro,m,) (24) 

quark 

The quantity Eym{R) is divergent due to singularity of Ad at r = R . We 
assume that Eym{R) become finite after renormalization. However, so far we have 
not elaborated renormalization procedure in our model. Instead, we simply postulate 
that 

EyMiR) = BR'' (25) 

where B is some constant. The choice n = 3 correspond to MIT model. In our model 
the choice n = 2 also seems natural because the main contribution in Eym{R) is 
proportional to the area of this sphere. 

Fortunately, it appears that hadronic masses depend on n very weakly. So concrete 
choice of this parameter is not important. 

Results of calculation of hadronic masses, corresponding to ground state configu- 
rations of quarks, are presented in Tables 3 and 4. Parameters -B,m„ = m^,ms, 
and rrih are determined by minimizing of the quantity 

A{B, m., m„ m„ m,) = ^[(K''' " Mf-O/^r]' 

h 

where M^^^ are masses of hadrons |T5], Mfl''^°' are defined by eq.(24). 

To be exact, for determination of the constant B and quark masses we used masses 
of p{n), S, Ac, A°, Df, Df^, B and rjc that are measured with the best accuracy in 
comparison with ones of other hadrons. 



5 Discussion. 

We see that our model gives rather good description for almost all hadron masses 
corresponding to ground states configurations of quarks except those of light pseu- 
doscalar mesons. As we already have mentioned in Introduction, the description of 
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Parameters 


n = 3 


n = 2 


B 


3.338 X lO^MeV)'^ 


2.960 X lO^MeVf 


m„ 


89.92{MeV) 


95.56{MeV) 


rris 


359.1{MeV) 


350A{MeV) 


rUc 


U32{MeV) 


U09{MeV) 


rub 


4893(Mey) 


m6(MeV) 



Table 3: Parameters of our model. 

the light pseudoscalar mesons cannot be given without consideration of the problem 
of chiral symmetry breaking that is out of the scope of our model nowadays. 

The accuracy 3-7 percents achieved in our model is maximal possible one for 
any constituent quark model in which interaction between quarks is not taken into 
account. Indeed, in any such model A and S particle must have the same mass. But 
really there exists approximately 7% difference between m/^ and m^.The same is true 
for Ac and S^. So seven percents is, most likely, the maximal accuracy that can be 
achieved in any simple constituent quark model. In fact, this means that we cannot 
describe spectrum of hadron resonances in framework of our model now. Indeed, 
mass differences between hadron resonances arc less, typically, then seven percents. 

Our model meets some internal difficulties. In particular, we are obliged to impose 
"by hands" boundary condition (8) for quarks with J > 3/2. May be this difficulties 
(as many others) can be overcome if one consider instead of the solution (6) more 
general solution of SU(3) YM equation. Indeed, in the present paper we investigate, 
in fact, SU(2) QCD with quarks in adjoint representation. But qauge group of QCD 
is definitely SU(3) rather than SU(2). Most likely, many difficulties arising in our 
model are connected with this circumstance. 

In the nearest future we plan to investigate singular solutions of SU(3) YM equa- 
tions that cannot be reduced to any solution of SU(2) ones and also to develop per- 
turbative theory on the background of such singular solutions. Results of the present 
paper show that there exists a good chance to obtain satisfactory description of mass 
spectrum and other properties of hadron in this way. 
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lO-u -L»J(i ^Jb 


^601 




Sc(J= 1/2, m = 2465) 


\Su ^Sg ll^c 


2326 


2318 


L)±(j = 0,m = 1869) 




1874 


1883 


(J = 0,m = 1968) 


l^s - l^c 


2029 


2033 


S(J = 0,m = 5278) 




5313 


5320 


r^c (J = 0,m = 2979) 


l-5e - IS, 


3011 


3000 


p(j= l,m = 769) 


\Su ^ 1>S'„ 


701 


721 


is:* (J = l,m = 892) 


I'S'^j ^Ss 


869 


887 


0(J= l,m = 1019) 


ISs - ISs 


1032 


1047 


T(J= l,m = 9460) 


ISb - ISb 


9856 


9830 


A++(J = 3/2,m = 1232) 


ISu ~ ^Su — iPu 


1191 


1159 


fi-(J = 3/2,m = 1672) 


ISs - IS, - IPs 


1671 


1646 


O0(J = l/2,m = 2700(?)) 


ISg — ISg — IS'c 


2491 


2480 



Table 4: Hadron mass spectrum evaluated by formulas (23)- (25). ISu — — \Ss 
means configuration in which two u quarks and one s quark in ground states with 
J = 1/2, etc; \Pu and IP^ means ground states of u and s quark with J = ?>/2. One 
notes that configurations l^Su — IS'u — l^S^ and ISg — lSs — lSs are forbidden by Pauli 
principle. So just configurations ISu — ^Su — iPu and ISs — ^Sg — IPs are true ground 
state configurations for A++ and The value of spin in these cases automatically 
equal to 3/2. 
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